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ABSTRACT 

An  experimental  study  of  frequency  spreading  in  acoustic  signals 
forward  scattered  from  model  tank  wind  driven  water  surfaces  has  been  per- 
formed. The  study  employed  wide  band  probing  signals  to  doubly  sample  the 
surface  scatter  process.  A high  degree  of  precision  and  demonstrated  repeat- 
ability was  produced.  Parameters  varied  In  the  study  included  grazing  angle, 
surface  roughness,  wind  direction,  and  acoustic  wavelength.  Rayleigh  rough- 
ness parameters  varying  from  0.7  to  12.  were  studied.  The  results  of  the 
study  consist  of  a large  number  of  measured  frequency  spreading  functions  and 
an  analysis  of  the  frequency  shift  of  the  mean  scattered  power  and  RHS  spread 
as  a function  of  wavelength  for  each  condition  of  wind  and  grazing  angle. 

This  analysis  Indicates  qualitative  agreement  with  current  theory  in  the  gen- 
eral form  of  the  spread  spectrum  at  high  and  low  roughness,  and  in  the  linear 
dependence  of  the  mean  shift  on  frequency.  The  data  showed  a marked  sideband 
asymmetry  in  symmetric  geometric  configurations,  however,  which  Is  not  consis- 
tent with  any  published  analytical  model.  In  addition,  at  short  wavelengths, 
frequency  spread  seemed  not  to  increase  linearly  with  carrier  frequency. 


J 


Introduction: 

In  recent  years  a number  of  studies,  both  theoretical  and  empirical, 
have  attempted  to  explain  the  nature  of  the  frequency  spreading  undergone  by 
acoustic  signals  as  a result  of  interaction  with  the  ocean  surface.  The  prob- 
lem is  of  considerable  Interest  in  the  design  of  long  range  sonar  systems  and 
in  the  estimation  of  sea  surface  statistics  by  acoustic  methods.  Recent  the- 
oretical work  has  Included  approximate  Integral  formulations  by  Roderick  and 

Cron^,  Parkins^,  Williams^,  Scharf  and  Swarts^,  and  Tuteur^,  a perturbation 
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theory  based  model  of  Harper  and  LaBlanca  ’ , and  an  analysis  by  Middleton 
based  on  a model  of  the  surface  as  an  ensemble  of  scatterers. 

Each  of  these  models  is  of  course  based  on  somewhat  different  assump- 
tions and  approximations.  It  is  not  surprising,  therefore,  that  the  predl  - 
tlons  made  by  them  are  not  in  complete  agreement.  There  seems  to  be  wide  ac- 
ceptance of  the  assertion  that  in  general  the  temporal  spectrum  of  wavehelghts 
in  the  area  of  the  surface  involved  in  the  scattering  will  be  a primary  factor 
in  the  determination  of  the  spreading  function.  There  is  no  agreement,  however, 
on  the  nature  and  magnitude  of  asymmetry  between  upper  and  lower  sidebands  to 
be  expected  and  on  the  question  of  the  appearance  of  second  and  higher  order 
harmonics  of  the  surface  height  spectrum  in  the  frequency  spreading  function. 

In  addition,  the  relative  Importance  for  frequency  spreading  of  diffraction 
effects  and  the  sensitivity  of  the  function  to  implicit  or  explicit  assump- 
tions about  the  asymmetry  of  surface  slopes  remains  an  open  question. 

Fortunately,  two  recent  experimental  studies  performed  in  the  ocean 

have  been  able  to  shed  some  light  on  these  questions  for  long  and  short  range 

9 

situations.  Brown  and  Frisk  performed  a number  of  measurements  in  Che  Pacific 
at  ranges  of  488  and  152  meters  (1600  and  500  ft.)  and  at  relatively  low  acous- 
tic frequencies  (100  to  500  Hz.)  and  wavehelghts  (29  cm.  RMS).  These  measurements 
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represent  single  bounce  weak  scattering  situations.  The  results  of  these  mea- 
surements were  somewhat  inconclusive.  Evidence  was  presented,  however,  of 
sideband  asymmetry  and  the  presence  of  possible  higher  order  functions  of  the 
surface  spectrum  in  the  acoustic  spectra.  The  asymmetry  was  confined  to  rela- 
tive heights  of  individual  components  and  total  sideband  power  was  in  general 
found  to  be  aymmetric. 

Shqoter  and  Mltchell^^  in  subsequent  measurements  at  much  longer 
ranges  (93  to  370  km.)  in  the  deep  ocean  and  at  substantially  lower  frequencies 
(92  and  137  Hz.)  gave  an  idea  of  the  spreading  to  be  expected  in  multiple 
bounce  RSR  situations.  Unfortunately,  this  aspect  of  the  experimental  situa- 
tion restricts  the  extent  to  which  the  results  can  be  used  as  a benchmark  for 
theory.  However,  the  presence  of  higher  order  structure  in  the  sidebands  and 
a marked  asymmetry  in  power  favoring  the  higher  sideband  observed  at  higher 
sea  states  (2.1  m.  sea  and  swell,  confused  surface)  reinforces  the  notion  that 
a simple  function  of  the  surface  spectrum  is  probably  not  sufficient  to  ex- 
plain frequency  spreading  due  to  surface  scatter. 

The  situation,  then,  is  that  several  mathematical  models  are  avail- 
able to  predict  frequency  spreading,  each  exhibiting  different  properties,  and 
no  experimental  evidence  is  available  with  which  to  sort  them  out.  The  pur- 
pose of  the  work  presented  in  this  paper  was  to  produce  experimental  data  of 
sufficient  precision  and  variety  to  permit  a more  decisive  relative  evaluation 
of  the  current  models.  Although  this  work  does  not  unilaterally  support  any 
of  the  models  mentioned  above,  it  does  clearly  show  a number  of  Che  controver- 
sial effects  in  relatively  clearly  defined  situations. 


The  Experimental  Method; 


The  experiments  discussed  in  this  paper  were  conducted  in  a model 
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tank  7.3  m.  In  diameter  and  1 m.  deep.  The  facility  used  has  been  described 
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previously  ' and  remains  largely  unchanged.  The  rough  surfaces  were  pro- 
duced by  a fan  driven  wind  tunnel  with  an  effective  fetch  of  3.4  m.  upwind  of 
the  scattering  area  and  wind  speeds  of  5.4  and  8.3  m./sec.  Reflections  of 
surface  waves  from  the  downwind  end  of  the  tank  are  damped  by  a "beach"  so  as 
to  suppress  free  gravity  waves  (swell) . This  was  done  so  as  to  study  the  ef- 
fects of  wind  driven  waves  alone.  All  measurements  were  made  with  projector  and 
hydrophone  directed  at  the  mean  surface  specular  point  and  at  equal  depth  and 
grazing  angle.  The  geometry  Is  defined  in  Figure  1. 

A technique  of  doubly  sampling  the  response  of  the  surface  to  a 
broadband  pulse  was  employed  In  this  series  of  measurements  as  a way  of  ac- 
quiring a large  amount  of  precise  data  efficiently  while  suppressing  spurious 
reflections.  The  procedure  is  based  on  the  preemphasized  pulsing  technique 
described  In  [11].  A preemphasized  probing  signal  s(t)  Is  designed  and  syn- 
thesized. This  signal  is  applied  to  the  projector.  Interacts  with  the  surface. 
Is  received  by  the  hydrophone  and  digitized  at  the  same  sampling  rate  used  to 
synthesize  it.  Because  of  the  linearity  of  the  process.  It  is  assumed  and  has 
been  demonstrated  that  the  signal  waveform  received  after  digitization  is  the 
convolution  of  the  original  probing  signal,  the  Impulse  responses  of  the  two 
transducers,  and  the  time  varying  Impulse  response  of  the  surface.  Since  the 
surface  Impulse  response  Is  slowly  varying  with  respect  to  the  time  scale  of 
the  signal  Itself,  It  Is  possible  to  represent  this  process  after  transforma- 
tion to  signal  frequency  domain  as  follows: 

R(u)^,t)  - S(u^)PQ(w)H(u,t)Pj^(«)  (1) 

where  Is  the  probing  signal  applied  to  the  projector 
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256  samples 


Pq(w)  Is  the  projector  transfer  function 

Pj^(u))  is  the  receiver  transfer  function 

H(u>,t)  is  the  time  varying  T.  F.  of  the  surface 

K(u  ,t)  is  the  signal  out  of  the  hydrophone  and  amplifiers 

13 

The  data  collection  instrumentation  removed  a fixed  time  delay 
from  the  received  signals  which  is  reflected  in  the  notation  here  in  the  use 
of  the  same  t in  both  H(a),t)  and  R((i)^,t)  and  the  suppression  of  time  depen- 
dence in  S(ui  ).  The  representation  of  S(a)  ) and  R(u)  ,t)  as  defined  on  a finite 
n n n 

domain  of  frequencies  is  based  on  the  fact  that  both  waveforms  are  wholly  con- 
tained in  relatively  short  time  windows  and  are  analysed  in  terms  of  the  dis- 
crete Fourier  transform. 

In  this  representation,  H(a),t)  is  a stationary  ergodlc  random  pro- 
cess expressing  the  amplitude  and  phase  response  of  the  surface  to  frequencies 
b)  at  time  t.  Care  was  taken,  of  course,  to  sample  in  t at  rates  in  excess 
of  the  Nyqulst  rate  of  the  process;  the  actual  sampling  rate  used  varied  be- 
tween 50  and  300  per  second.  In  practice,  sampling  in  t was  limited  to  64 
contiguous  samples  in  a burst.  This  limit  was  imposed  by  memory  availability 
in  the  data  acquisition  computer,  and  is  equivalent  to  CW  pulse  lengths  of 
from  850  to  20,000  cycles  of  the  carrier.  The  recorded  samples  of  H(u,t)  were 
then  transformed  on  the  variable  t to  give  an  ensemble  of  Instantaneous  spread- 
ing functions. 

^ N-1  -Jv.mAt 

■'‘"n-V  ■ ) (2) 

m-0 


An  ensemble  of  roughly  200  measurements  of  were  then  windowed  and 

averaged  to  give  the  frequency  spreading  function. 
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where  frequency  domain  Hanning  window.  The  design  of  the  signal 

S(u^^)  was  such  Chat  Che  first  term  in  (2)  was  approximately  unity.  This  re- 
sulted in  a uniformly  precise  estimate  of  over  the  frequency  domain 

of  Che  probing  signal. 

This  procedure  turned  out  to  have  substantial  advantages  over  the 
alternative  of  using  single  long  CU  pulses.  The  advantages  are  the  near  total 
immunity  from  multipath  interference  made  possible  by  time  gating  of  broad 
band  pulses,  and  parallel  acquisition  of  data  at  a large  number  of  carrier 
frequencies  simultaneously.  The  principal  disadvantage  at  present  is  the  rel- 
atively restricted  resolution  in  Doppler  frequency  obtainable  with  the  present 
available  memory.  This  turned  out  not  to  be  a problem  in  this  series  of  mea- 
surements . 

Data  were  collected  using  this  scheme  for  two  wind  speeds,  two  graz- 
ing angles,  and  two  orientations  with  respect  to  Che  wind  direction.  In  each 

case,  over  100  values  of  u were  available  from  Che  data  of  which  4 values 

n 

were  selected  for  presentation.  These  results  are  displayed  in  Appendix  A as 
plots  of  Doppler  frequency  vs.  relative  power  for  a variety  of  values  of 
wind  speed  and  geometry.  The  water  surface  temporal  spectra,  measured  by  a 
capacitance  probe,  for  each  of  the  two  wind  speeds  are  shown  in  Figs.  3 and  4. 
The  general  properties  of  Che  short  fetch  wind  driven  water  surfaces  produced 
in  the  model  tank  have  been  discussed  previously^^'^^'^^.  The  spe'tra  have 
been  included  here  because  of  their  relevance  to  several  of  Che  theories  re- 
ferred to.  These  measurements  represent  conditions  of  strong  scattering, 
crosswind  and  downwind  orientation,  and  long  and  short  wavelengths  with  respect 
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Fig.  3: 

Surface  tciaporal  power 
spectrum  for  Wind  2, 
o * 1.1  an. 
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to  the  surface  roughness.  In  all  cases,  the  beam  patterns  of  the  projector/ 
receiver  pair  were  substantially  larger  than  the  first  several  Fresnel  zones 
on  the  surface. 

As  a check  on  the  linearity  and  repeatability  of  the  data  acquisi- 
tion apparatus,  an  experiment  was  performed  which  was  In  all  ways  Identical 
to  that  of  Fig.  A-1,  except  for  a modification  of  the  probing  signal.  The 
signal  used  to  obtain  the  spectra  in  Appendix  A was  designed  to  present  the 
receiver  with  a flat  spectrum  pulse  having  essentially  equal  energy  at  all 
frequencies  between  0.2  and  1.2  MHz  (Fig.  5).  The  pulse  used  in  the  test  run 
was  similar  except  for  deletion  of  spectral  components  between  0.6  and  0.9 
mHz  (Fig.  6).  These  two  signals  obviously  differ  greatly  in  their  time  do- 
main forms.  Therefore,  they  could  be  expected  to  provide  a good  test  of  the 
linearity  assumptions  upon  which  the  data  analysis  procedure  used  depends. 

Fig.  7 shows  the  results  of  the  test  measurement  for  256  kHz.  and  1.06  MHz. 

As  can  be  seen  by  comparison  with  Fig.  A-1,  the  degree  of  agreement  Is  quite 
good. 

Discussion; 

Inspection  of  the  results  presented  In  the  previous  section  reveal 
the  following  qualitative  observations: 

1.  Asymmetry  of  frequency  spread  is  observed  to  be  the  rule  rather 
than  the  exception,  even  In  the  symmetric  geometric  configuration.  The  amount 
and  direction  of  the  asymmetry  appears  to  be  dependent  on  all  of  the  parameters 
varied  In  the  study.  At  low  Rayleigh  parameters,  sidebands  appear  to  have  un- 
equal power  although  the  locations  of  peaks  in  the  sideband  spectra  are  coin- 
cident with  peaks  In  the  wavehelght  temporal  spectra.  This  phenomenon  Is 
wholly  consistent  with  the  notion  of  a system  function  as  used  by  Harper  and 
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LaBianca  , Williams  , and  others.  Its  presence  in  a symmetric  geometry,  how- 
ever, indicates  that  there  are  significant  inadequacies  in  the  formulations 
presently  used  to  calculate  these  functions.  The  transition  to  high  roughness 
appears  to  take  place  between  Rayleigh  numbers  1 and  2.  In  this  domain,  a 
relatively  smooth  approximately  Gaussian  spectrum  appears  to  add  to  the  weighted 
image  of  the  surface  spectrum  until  at  values  greater  than  2 it  dominates  the 

frequency  spreading.  In  this  regime,  the  nature  of  the  spread  spectrum  appears 

4 

to  be  generally  consistent  with  the  prediction  of  Scharf  and  Swarts  in  that 
the  spread  spectrum  is  generally  Gaussian  and  with  a mean  shifted  from  the 
carrier  by  an  amount  proportional  to  the  carrier  wave  number.  However,  this 
model  also  predicts  symmetry  in  the  configurations  used  here,  and  therefore 
does  not  completely  account  for  the  observed  spectra.  A suggested  reason  for 
this  disagreement  will  be  given  later  in  this  paper.  The  linear  dependence 
on  wave  number,  however,  is  consistent  with  the  data,  at  least  qualitatively 
(Fig.  8). 

2.  The  frequency  spreading  bandwidth  increases  with  Increasing  car- 
rier frequency,  but  not  linearly  end  not  without  limit.  RMS  spreading  band- 
width versus  carrier  frequency  is  plotted  in  Figure  9.  In  the  regime  of  low 
roughness,  the  bandwidth  seems  to  Increase  slowly  from  an  asymptotic  value 
equivalent  to  the  bandwidth  of  the  surface  temporal  spectrum.  This  is  hardly 
surprising  and  consistent  with  the  models  valid  in  that  domain.  In  addition, 
the  data  of  Roderick  and  Cron,  Brown  and  Frisk,  and  Williams  reflect  a similar 
behavior.  Roderick  and  Cron  have  proposed  an  explanation  in  terms  of  higher 
order  sidebands  due  to  a simultaneous  phase  and  amplitude  modulation  process. 

At  higher  values  of  roughness,  however,  Scharf  and  Swarts  predict  a linear 
dependence  on  carrier  frequency  which  is  not  supported  by  the  data.  Indeed, 


race  of  Increase  of  bandwidth  with  carrier  frequency  appears  to  caper  off  Co 
zero  and  in  some  cases  decreases.  Apparently,  Che  point  at  which  Che  spread- 
ing stops  increasing  is  not  a function  of  surface  roughness  or  wind  direc- 
tion, but  only  of  grazing  angle  and  wavelength.  This  has  led  the  author 
to  speculate  that  the  cause  of  the  levelling  off  may  have  to  do  with  the  de- 
crease in  Che  size  of  the  first  several  Fresnel  zones  on  Che  surface  which 
are  Che  principal  contributors  to  the  signal  received  from  the  surface. 

3.  Spreading  in  Che  crosswind  orientation  is  substantially  stronger 
than  in  Che  upwind  direction.  To  a qualitative  extent,  this  effect  is  consis- 
tent with  the  results  of  Scharf  and  Swarts  Chat  the  spread  can  be  predicted 
from  Che  second  derivatives  of  the  surface  space-time  correlation  function  at 
the  origin.  Quantitative  verification  would  require  both  an  extensive  set  of 
measurements  of  the  surface  correlation  functions  at  small  lags  and  displace- 
ments, and  a resolution  of  the  symmetric  geometry  issue  referred  to  above. 

The  author  is  considering  some  practical  approaches  to  both  problems. 

Middleton  has  suggested  that  the  Inadequacy  of  current  treatments 
of  the  Doppler  spreading  process  stems  from  the  tendency  of  the  various  treat- 
ments to  incorporate  incomplete  treatmenCs  of  the  two  basic  mechanisms  of  fre- 
quency spreading.  Doppler  spread  in  signals  reradlated  from  a rough  moving 
surface  arises  both  from  individual  physical  interactions  of  the  incident 
field  with  moving  elements  of  the  surface  (phase  modulation),  and  from  the 
interference  of  contributions  to  the  reradlated  field  due  to  variations  in 
path  length  (amplitude  and  phase  modulation).  Surely  these  two  mechanisms  are 
intimately  coupled  through  the  dynamics  of  the  wind  driven  water  surface. 
Nevertheless,  in  all  of  the  analytical  treatments  cited  above,  simplifying 
assumptions  have  been  made  which  have  the  net  effect  of  either  inaccurately 


modeling  one  mechanism  or  the  other,  or  of  neglecting  the  hydrodynamic  coupling 


between  the  two.  In  the  data  presented  In  this  paper,  as  well  as,  to  a lim- 


ited extent.  In  the  ocean  data,  a substantial  difference  can  be  seen  between 


high  and  low  roughness  spreading.  This  difference  Is  not  only  In  magnitude 


but  also  In  form  and  may  well  have  to  do  with  a change  In  the  relative  Impor- 


tance of  the  two  effects. 


Another  obvious  and  universally  recognized  weakness  of  the  analyt- 


ical models  Is  In  the  relatively  simple  expressions  usually  used  for  the  sea 


surface.  Harper  and  LaBlanca  use  a fairly  typical  expansion  In  sinusoids  of 


the  form: 


Z(^,t)  - ^ I E h^jCos(r^j  • p - fijt  + 


where 


Irijl  -lijg 


and  the  component  phases,  assumed  to  be  uniformly  distributed  and 


statistically  Independent.  This  assumption  of  uniform  Independent  phases  Is 


probably  not  realistic  but  is  mathematically  very  convenient.  Scharf  and 


Swarts  use  an  equivalent  representation  expressed  In  terms  of  a spectral  den- 


sity. They  correctly  point  out  that  this  representation  is  valid  only  for  a 


restricted  class  of  Gaussian  surfaces  which  are  not  significantly  wind  driven 


or  turbulent.  In  view  of  the  relatively  well  understood  nature  of  scattering 


from  sinusoidal  Interfaces,  such  an  expansion  has  obvious  advantages.  Indeed, 


for  a surface  dominated  by  swell,  or  for  a mechanically  excited  surface  such 


as  the  one  studied  by  Roderick  and  Cron  , this  form  Is  a reasonably  good  repre- 


sentation. However,  for  surfaces  in  which  a substantial  contribution  is  due 


to  the  local  wind  (sea) , the  sinusoidal  model  with  random  phases  Is  certainly 
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inaccurate.  The  most  convincing  evidence  of  this  is  to  be  found  in  measure- 
ments of  water  surface  slope  distributions^^’^^’^^  and  wind  direction  depen- 
dence of  backseat ter^^.  The  surfaces  produced  in  the  model  tank  used  in  this 
study  have  been  designed  to  have  a very  minimal  swell  and  are  predominantly 
wind  driven.  It  should  not  be  surprising,  therefore,  that  the  Doppler  spread- 
ing observed  is  not  entirely  predicted  by  the  analytical  models  which  have 
been  based  on  an  assumption  of  a simple  surface  representation.  Actual  wind 
driven  surfaces  are  very  difficult  to  model  in  a realistic  fashion.  However, 
in  order  to  produce  a generally  useful  predictive  theory  for  ocean  scattering, 
some  improvements  in  the  present  surface  representations  will  probably  be 
required . 


Conclusion; 

The  results  of  this  study  indicate  two  general  conclusions: 

1.  The  current  theories  proposed  for  the  prediction  of  frequency 
from  rough  water  surfaces  do  to  some  extent  work.  Th&y  predict  a modified 
image  of  the  surface  temporal  spectrum  in  the  case  of  small  roughness,  and  a 
shifted  Gaussian  of  greater  width  in  the  case  of  large  roughness.  Both  of 
these  predictions  seem  to  be  qualitatively  consistent  with  the  experimental 
measurements  on  wind  driven  surfaces. 

2.  In  general,  the  quantitative  predictions  of  the  theories  are 
not  accurate  fur  the  case  of  wind  driven  surfaces.  The  most  obvious  disagree- 
ments are  that  the  observed  spreading  is  not  symmetric  under  conditions  of 
symmetric  geometry,  and  that  the  bandwidth  of  the  scattered  sound  does  not 
increase  linearly  with  acoustic  frequency.  The  exact  reasons  for  this  dis- 
agreement are  at  present  obscure,  but  seem  likely  to  be  the  result  of  Incomplete 
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surface  representations  In  the  analyses.  These  conclusions  Indicate  that  the 
theoretical  models  are  In  general  correctly  formed,  but  require  additional 
refinement  In  order  to  be  useful.  The  author  Intends  to  direct  his  attention 
to  measurements  on  surfaces  with  substantial  swell,  off-specular  scattering, 
and  studies  of  the  second  order  statistics  of  the  surfaces.  In  this  way.  It 

; 

may  be  possible  to  arrive  at  a useful  compromise  between  the  limited  generality 
of  experimentation  and  the  necessary  evils  of  analytical  approximation.  | 
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Appendix  A; 


The  following  plots  are  of  the  measured  frequency  spreading  for  a 
variety  of  conditions  and  carrier  frequencies.  Each  is  identified  by  an  iden- 
tification block.  The  information  in  the  block  is  given  in  the  following 
format : 

Run  No.  Wind  No. 

XMT  RQ(cm)  iJ/^Cdeg)  8Q(deg) 

REC  Rj^(cm)  i))^(deg)  6^(deg) 

Carrier  Freq. 

Surface  roughness  is  commonly  expressed  as  either  RMS  wavehelght  or  Rayleigh 
delay.  These  values  are: 


Wind  No. 

a (mm. ) 

4ira  sin  , s 

(sec) 

c 

2 

1 

o 

1.13 

2.8  X lO"^ 

2 

30° 

1.13 

4.8  X 10“^ 

4 

17° 

2.6 

6.5  X 10"^ 

4 

0 

O 

2.6 

11.1  X 10“^ 
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Figure  A-4a 
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Figure  A- 5b 
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Figure  A-8b. 
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